Mon. Not. R. Astron. Soc. 000,[I]-?? (2010) Printed 15 November 2010 (MN WT^ style file v2.2) 



A Wide Field Survey of Satellite Galaxies around the 
Spiral Galaxy M106 
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ABSTRACT 

We present a wide field survey of satellite galaxies in M106 (NGC 4258) 
covering a 1.7° x 2° field around Ml 06 using Canada- France-Hawaii Tele- 
scope/MegaCam. We find 16 satellite galaxy candidates of M106. Eight of 
these galaxies are found to be dwarf galaxies that are much smaller and 
fainter than the remaining galaxies. Four of these galaxies are new find- 
ings. Surface brightness profiles of 15 out of 16 satellite galaxies can be 
represented well by an exponential disk profile with varying scale length. 
We derive the surface number density distribution of these satellite galax- 
ies. The central number density profile (d < 100 kpc) is well fitted by a 
power-law with a power index of —2.1 ± 0.5, similar to the expected power 
index of isothermal distribution. The luminosity function of these satel- 
lites is represented well by the Schechter function with a faint end slope 
of — l.lO^Q Qg. Integrated photometric properties (total luminosity, total 
colour, and disk scale length) and the spatial distribution of these satellite 
galaxies are found to be roughly similar to those of the Milky Way and 
M31. 
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1 INTRODUCTION 

Dwarf galaxies are believed to be the most abun- 
dant population in the universe. The preponderance 
of dwarf galaxies has been pr edicted by co l d dar k 
matter (CDM) simulations feg.. lKlvDin et al.1 (|l999l l) 
but the slope of the observed luminosity function of 
faint galaxies is too shallow to be matched with the 
predictions of CDM cosmology. Another important 
prediction of hierarchical structure formation in the 
CDM or ACDM cosmology is that most of the dwarf 
galaxies are thought to be satellites of giant galax- 
ies because galaxies are as sembled b y sub-structures 
of sub-galactic masses (^Frenk et al.1 11988 ) and many 
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low-mass haloes can survive after major and minor 
mergers. 

To date there are only a few giant spiral galax- 
ies which are known to have more than dozens 
of satellite galaxies: the Milky Way Galaxy (MW), 
M31 and the M81 group. The number of satel- 
lite galaxies in MW and M31 has been dou- 
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120071 ). mostly from the Sloan Digital Sky Survey 
(SDSS) (lYork et al.l l2000V Three of them (Willman 
1, Segue 1, Bootes II) are very faint [My ^ — 3) and 
it is uncertain whether they are genuine dwarf galax- 
ies ('Willman et all l2005bl : iBelokurov e\~^ l2007bl : 
■Walsh et al. 2007 ) . The number of satellite galaxies in 



© 2010 RAS 



2 E. Kim et al 



Table 1. Basic information of M106 (NGC 4258) 



Parameter Values References 



RA(2000), Dec(2000) 12hl8m57.5046s, +47° 18' 14'^303 1 

Ellipticity 0.611 2 

RA. 150 deg 2 

Standard radius, R25 9.3 arcmin 2 

Total magnitudes = 8.41 ± 0.08, = 9.10 ± 0.07 2 

X-ray luminosity Log(Lx) = 40.88 erg s"-*^ 3 

Heliocentric radial velocity, Vr 448 zb 3 km s"-*^ 1 

Foreground reddening E{B - V) = 0.016 4 

Distance d=7.2 Mpc ((m - M)o = 29.29 ± 0.07) 5,6 



(1) NASA Extragalactic Database (NED); (2) RC3; (3) Panessa et al (2006) (nuclear 
flux - Chand r a obs ervation); (4) Schlegel et al. (1998); (5) Herrnstein et al. (1999); (6) 
iMager et aP (l2Q08l ) 



the M81 group has been doubled too by the recent dis- 
cover ies of 22 new satellite galaxies ([Chiboucas et all 
I2OO9I ) from a deep survey of 8° x 8° area using Mega- 
Cam of CFHT. Among them, the satellite galaxies of 
the M81 group are of particular interest because they 
are faint enough to be compared with those of MW 
and M31, and provide unique information about the 
formation of satellite systems in a group environment 
that is thought to be different from that of the Lo- 
cal Group. Because of the very wide searching area of 
8° X 8° , correspon ding to 500 kpcx 500 kp c at the dis- 
tance of 3.6 Mpc, Chiboucas et al. ('2009) fairly well 
covered all the satellite galaxies within the virial ra- 
dius of M81. The dramatic increase of the number of 
satellite galaxies in these systems seems to greatly al- 
leviate the missing satellite problem. 

However, there is still a discrepancy between 
the number of observed dwarf satellite galaxies and 
those predicted by the CDM hierarchical models 
(|Klvpin et al.lll99d : iMoore et all 11999). because the 
number of observed MW satellites corrected for the 
sky coverage of the SDSS is st ill a factor of 4 smalle r 
than the model predictions (Simon k Gehal I2OO7I ). 
There have been several explanations to solve the dis- 
crepancy between the predicted and observed num- 
bers of dwarf galaxies. Some solutions modify the pri- 
mord i al density fluctuations (p^mionkowski Liddld 
I2OOQI : IZentner Bullockl |20Q3') or the dark matter 
prop e rties by invoking w arm dark matter (Coli n et al.l 
I2OOQI : iBode et al.1 l200lh . However, more appealing 
explanations employ astrophysical solutions to sup- 
press gas accretion on the low mass dark haloes or 
to prevent star formation in these haloes so that 
they can not be observed as dwarf galaxies. Among 
others photoionizati o n squelching afte r re-ionization 
iBullock et al. 1 I2OOQI : ISomervillelT2002l : iBenson et all 
[2002) seems to be the most att ractive because it pro - 
vides a natural low mass cutoff (Koposov et al.l l2009l ). 
However, iKravtsov et al. 1 C2OO4) argued that the star 
formation history in the low mass dwarf galaxies is 
not sensitive to UV background and re-ionization but 
sensitive to the physics of galaxy formation. They pro- 



posed tidal stripping as an alternative mechanism for 
the paucity of observable dwarf galaxies at the present 
epoch. 

Owing to the recent discoveries of new satellite 
galaxies in MW, M31, and the M81 group and some 
models that can explain the observed properties of 
satellite galaxies belonging to these galaxy haloes, the 
missing satellite problem seems to be less severe than 
before. However, satellite galaxies in the Local Group 
(LG) may not be a representative population of dwarf 
galaxies surrounding dark haloes because LG is lo- 
cated in a underdense region. Since the number of 
bright members of the M81 group is similar to that 
of LG and their local background densities seem to 
be not much different due to the very close distance 
between LG and M81. Thus, we need more sample in 
different environment to understand this problem. 

The purpose of the present study is to find dwarf 
satellite galaxies surrounding a giant spiral galaxy 
which is similar to MW and M31 to address the miss- 
ing satellite problem in different group environment. 
We selected M106 (NGC 4258) as the best target for 
our deep and wide imaging survey to search for satel- 
lite galaxies because of its proximity (^ 7.2Mpc) and 
brightness {My = —20.93) that allow us to detect 
faint satellite galaxies down to My ~ —10. Ml 06 
is the brightest member of CVn H group which is 
thought to be one of the most dense grou ps in the local 
univer s e according to the catalogue of iFouque et al.l 
("199211. iTurner GottI ([197 6) assigned M106 to their 
Group 50 of which some galaxies are members of CVn 
H. 

The morphology of Ml 06 is classified as a barred 
spiral galaxy (SAB(s)bc). M106 is known as a LINER 
harboring nuclear water masers which can be used 
for an accurate distance estimate by geometrical 
means. The distan ce estimated from the w ater maser 
is 7.2 zb 0.3Mpc ("Herrnstein e t al.1 Il999l l. which is 
consistent with the distances determined using the 
Period— Luniinosity relation of Cepheids and TRGB 
(^Macri e t af] I2OO6I: [Benedict et al.ll2007l : iMager et all 
■2008 : Bono et al.l I2OO8I I . At this distance one arcsec 
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Figure 1. A finding chart for the four fields of M106 marked on the Digitized Sky Survey image. The size of field of view 
is 2.2° X 2.3°. North is up and east is to the left. Satellite galaxies are marked by circles with running ids and the extent 
(D25) of M106 is shown by a large ellipse. Probable satellite galaxies PI— P5 of M106 are marked with open squares. 
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corresponds to a linear scale of 35pc. The luminosity 
of M106 is slightly fainter than those of MW and M31. 
Its systemic velocity is 448 ± 3 km s~^. Its maximum 
rotation velo city, 208 km s~\ is slightly smaller than 
that of MW (|Erickson et al.ll 19991 ). It is known to have 
one satellite galaxy, NGC 4248 (~ 13 arcmin from 
the centre of M106). Foreground reddening toward 
M106 is low: Ay = 0.05, Ar 0.04 and Ai 0.03 
fechl egel et al.lll998l l. Basic information on M106 is 
summarized in Table 1. 

We surveyed an area of 1.7° x 2° centered on 
M106, which corresponds to 215 kpcx 260 kpc, using 
g and r filters. Our surveyed area is expected to cover 
about half of the virial radius where most of the dwarf 
sa tellite galaxies are considered to reside. The model 
oflMoore et al.' (l999') predicted that about 75% of the 
dwarf satellite galaxies are located within r = 100 kpc 
from the centre of the host dark halo. Thus, we expect 
to detect more than 75% of the satellite galaxies from 
the present survey. 

In this paper we report on the candidate satellite 
galaxies of Ml 06 some of which are newly discovered 
dwarf galaxies from a deep MegaCam survey. Obser- 
vations and data reduction are described in §2. We 
present the main results of the present study in §3. 
In §4 we discuss the misdentifications found in the 
present study and missing satellite galaxies, and com- 
pare the satellite system of Ml 06 to those of LG and 
the M81 group. We summarize the main results of this 
study in the last section. 



2 OBSERVATIONS AND DATA 
REDUCTION 

2.1 Observations 

CCD images of M106 were obtained on the 4 nights 
from February 11 2005 to April 6 2005 (UT) using the 
wide-field camera, MegaPrime/MegaCam at CFHT. 
Most of observations were carried out under a good 
photometric condition except for the observations of 
March 17 2005 when it was cloudy. MegaCam is made 
of 36 2048 X 4612 pixel CCDs (4 rows x 9 columns of 
CCDs), and it covers 1° x 1° field of view with a reso- 
lution of 0.185 arcsec per pixel at CFHT MegaPrime 
focus. SDSS-^ and r filters were used for our observa- 
tion. Four fields around M106 (Fl, F2, F3, and F4) 
were observed and the total field coverage is 1.7^^ x 2^, 
as displayed in Fig. 1. All fields were observed to be 
partially overlapped with each other to check internal 
photometric accuracy. The observation log is given in 
Table 2. Exposure times are 4 x 210 s for g filter, 
and 4 x 600 s for r filter. Each image was taken with 
dithering of 20 ~ 30 arcsec. These dithered observa- 
tions turned out to be very useful to fill the gaps in the 
mosaic CCDs in the same row. The dithering pattern 
used in the present observation, however, could not 
cover the large gaps 0.9 arcmin) between different 
rows of CCDs, a gap between the first and the sec- 



ond rows and a gap between the third and the fourth 
rows. This uncovered sky region will reduce the ef- 
ficiency of finding satellite galaxies, especially small 
systems (smaller than 20^^). The seeing ranged from 
0.58 to 1.16 arcsec with a median value of 0.74 arcsec. 



2.2 Photometry 

The raw images were preprocessed using Elixir sys- 
tem by CFHT staff. Elixir is a collection of programs, 
databases, and other tools specialized in processing 
and evaluation of the large mosaic data. Detailed in- 
formation about the Elixir system can be found in 
iMagnier k Cuillandrd (|2004l ). 

We combined the dithered exposures of each CCD 
to get rid of bad pixels, defective columns and spuri- 
ous cosmic rays by taking a median of pixel values. 
We utilize the resulting images for the detection of 
point source object s with the dig ital photometry pro- 
gram DAOPHOT (|Stetson|[l994 ). We adopt 3a as a 
detection threshold. Instrumental magnitudes of the 
detected objects in the images are derived using the 
aperture photometry with aperture radii of 0.6 to 1.2 
arcsec depending on the seeing FWHM of a point 
source on each CCD. The value of the aperture cor- 
rection was derived from the difference between these 
small aperture magnitude and large aperture magni- 
tude for well isolated bright stars. 

To derive the radial surface brightness profiles 
of the satellite galaxies we perform surface photom- 
etry. We define the rough boundary of the satellite 
galaxies by estimating major axis radius, ellipticity 
and position angle for each satellite galaxy (4th, 5th 
and 6th columns in Table 3) by visual inspection. 
Then we obtain the radial surface brightness profiles 
of the satellite galaxies using a custom software El- 
lipsephot Ellipsephot performs a radial surface pho- 
tometry using elliptical apertures by dividing an el- 
liptical annulus into sectors to deal with unexpected 
spurious sources and/or bad columns. For the surface 
photometry of candidate satellite galaxies of Ml 06 we 
divide an elliptical annulus into 8 sectors and take 
the median value of surface brightness of these sectors 
as the surface brightness of certain radius. This pro- 
cess is useful to remove the effect of spurious sources 
and/or bad columns in the surface photometry of a 
galaxy. Preliminary version of Ellipsephot was success- 
fully applied in a study of the surfa ce photometry o f 
dwarf galaxies in the Local Group (|Lee et al.lll999l V 
Ellipsephot is useful to obtain the radial surface pho- 
tometry of non-elliptical shaped and/or very low sur- 
face brightness galaxies, because it is not involved with 
fitting process. Therefore, Ellipsephot can be used for 
the surface photometry of all extended sources ir- 
respective of their shapes. Users need to determine 
the centre, ellipticity, and position angle of the ex- 
tended structure to run Ellipsephot. For galaxies with 
elliptical structures we obtain their surface photom- 
etry using IRAF/STSDAS/ellipse task. IRAF/ellipse 
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Table 2. Observation log 



Field 


Filter 


Exp. Times 


Airmass 


Seeing^ 


Date(UT) 


Sky Brightness^ 


Lunar Angle^ 


Lunar Phase 


Fl 


g 


4 X 210s 


1.15 


0'^61 


2005 Feb 11 


22.42±0.05 


133 


0.15 


Fl 


r 


4 X 600s 


1.17 


0'^56 


2005 Feb 11 


21.58±0.03 


133 


0.15 


F2 


g 


4 X 210s 


1.19 


0'^62 


2005 Feb 16 


22.19±0.06 


95 


0.63 


F2 


r 


4 X 600s 


1.15 


0'^67 


2005 Feb 16 


21.50±0.03 


94 


0.63 


F3 


g 


4 X 210s 


1.18 


0'^78 


2005 Mar 17 


21.26±0.05 


79 


0.46 


F3 


r 


4 X 600s 


1.15 


0'^76 


2005 Mar 17 


20.94±0.02 


78 


0.46 


F4 


g 


4 X 210s 


1.15 


0'^95 


2005 Apr 06 


21.77±0.02 


140 


0.01 


F4 


r 


4 X 600s 


1.17 


o'^7l 


2005 Apr 06 


21.11±0.02 


140 


0.01 



"Typical seeing size at the centre of focal plane. Significant image quality degradation occurs in the outer 
regions due to the large-format CCDs. 
^In units of mag arcsec"^. 
•^In units of degree. 



task performs ellips e fitting to the observed isophotes 
(jjedrzeiewskill 19871 1 . A comparison of the surface pho- 
tometry for elliptical galaxy NGC 4346 (S15, see Table 
3) shows an excellent agreement between the results 
using Ellipsephot and that of IRAF/ellipse. 



2.3 Standard Calibration 

CFHT staff provided the transformation equations for 
a standard calibration from the photometry of the 
standard stars observed during the same night. We 
obtained the instrumental magnitudes of the standard 
st ars usin ^ the aperture radius of 7.5 arcsec as used 
m iGeisled (|l996l ). The standard transformation equa- 
tions are: g = g' + 0.148(^' - r) - 0.150X + const 
with rms=0.019 and N=72, and r = + 0.000(^' - 
r)-0.100X + const with rms=0.008 and N=66, where 
the lower case letters represent the standard magni- 
tudes, the primed lowercase letters the instrumental 
magnitudes (with a DAOPHOT system zero point of 
25.0), and X the air mass. Zero-points of these trans- 
formation relations are slightly different for different 
observed field. We transformed the instrumental mag- 
nitudes of the sources detected in fields Fl and F2 
onto the standard system using these transformation 
equations. The instrumental magnitudes of the point 
sources detected in F3 and F4 were transformed onto 
the standard system by comparing the magnitudes of 
the point sources common between Fl and F3 and 
between F2 and F4, respectively. 

SPSS Data Releas e 5 (DR5, 
lAdelman-McCarthv et all (|2007l )) covers the en- 
tire field of the present MegaCam observation of 
M106, though SDSS observation is significantly shal- 
lower than MegaCam observation. We find common 
point sources between SDSS DR5 catalogue and this 
study, and derive the magnitude difference between 
the transformed standard magnitudes of this study 
and the point-spread function fitting magnitudes of 
point sources in SDSS DR5. We adjust the standard 
magnitudes of the point sources by adding these 
magnitude differences in each CCD image. 



2.4 Sky Brightness 

Sky brightness is a critical parameter to find faint 
satellite galaxy. Detection efficiency of the extended 
sources (e.g., galaxies) strongly depends on both sur- 
face brightnesses of the sources and the background, 
while that of the of point sources depends mainly on 
the total source brightness. Therefore it is important 
to carry out observations on the dark night to find as 
many faint satellite galaxies as possible. 

We measure the sky brightnesses of the individ- 
ual fields, and list their values in Table 2. We select a 
few sky regions around the satellite galaxies to mea- 
sure the sky brightness located 3 — 5 times of the 
outer radii (see Table 3 for the values of the outer 
radii of the satellite galaxies) where no foreground 
and background sources are included. The mean sky 
brightnesses of Fl and F2 in g and r bands are 22.31 
and 21.54 mag arcsec"^, respectively, corresponding 
to 22.78 and 21.86 mag arcsec"^ in B and V bands. 
We transform the SDSS g and r magnitudes to John- 
son B a nd V magnitudes usin^ th e transformation re- 
lation in Chonis & Gaskell [2008) , derived by compar- 
ing SDSS DR5 magnitude of standard stars of Lan doltl 
(1992) to their standard magnitude. The derived sky 
brightness is slightly fainter than the typical dark sky 
brightness at the CFHT site (22.00 and 21.40 mag 
arcsec"^ in g and r bands). The mean sky brightness 
of F3 is significantly brighter than those of Fl and F2 
by 0.8 and 0.6 mag in g and r bands, respectively. Sky 
brightness of F4 is slightly fainter than that of F3, 
but is brighter than those of Fl and F2 by ^ 0.5 mag. 
The cause of the relatively brighter mean sky bright- 
ness of F3 compared to Fl and F2 might be due to the 
scattered moon light by the cloud and rather smaller 
lunar angle, and larger lunar phase. The extinction by 
cloud for F3 might affect the detection limiting mag- 
nitude and limiting surface brightness, therefore, the 
detection efficiency of satellite galaxies in F3 might be 
lower than those of the other fields. We list the lunar 
angle and the lunar phase of each observation field in 
Table 2. 
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^ [Sir, [52] NGC4218 [S3] NGC4226 : [$4^ NGC42il8 ^* 




Figure 2. An atlas of the satellite galaxies in M106. The numbers at the bottom of individual panels represent the size of 
the field of view. 



3 RESULTS 

3.1 Selection of Satellite Galaxies 

Distinguishing the satellite galaxies from the back- 
ground galaxies and foreground objects is one of essen- 
tials in the present study. Based on the simple hypoth- 
esis of similarity in physical properties of the satellite 
galaxies of Ml 06 and the Local Group, the satellite 
galaxies of Ml 06 are expected to have a variety of size, 
total magnitudes, colours and surface brightnesses. 
We therefore, need to explore a large set of searching 
parameters to complete the finding of satellite galaxies 
of Ml 06. Most efficient and complete method suitable 



for finding satellite galaxies in the wide field survey is 
human visual inspection. Three of the authors search 
both the combined fits images and true color images 
independently and repeat the searching process until 
three authors agree one another. The galaxies identi- 
fied as satellite galaxies by either one or two authors 
are classified as probable satellite galaxies. 

All the extended sources detected in the im- 
ages might be potential satellite galaxy candidates of 
Ml 06. However, considering the total mass (similar 
to or smaller than the mass of Milky Way), M106 
is expected to host a few tens of satellite galaxies. 
Therefore many of the extended sources in the im- 
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Figure 3. Surface brightness profiles of the candidate sateUite galaxies in M106 against the semi-major axis radii. The 
curved broken lines represent the point spread function. Locations of the effective radii, standard radii and Holmberg radii 
are marked with vertical arrows, vertical solid lines, and vertical dotted lines, respectively. The r-band surface magnitude 
is shifted by -0.5 mag for easier comparison. Surface brightness profiles are fitted by the exponential law except for S15, 
whose profiles are fitted with de Vacouleurs law. The curved solid lines over the data points represent the best-fitting with 
the exponential law or de Vacouleurs law. 



ages are not the satellite galaxies of Ml 06. To select 
genuine satellite galaxies while keeping as many can- 
didate galaxies as possible, we construct the following 
steps to find the satellite galaxies. 

(1) If a galaxy that shows some resolved stars, it is 
considered to be a satellite galaxy. Since we can resolve 
a significant portion of Ml 06 into individual sources, 
the partially resolved galaxies near Ml 06 could be the 
most plausible candidates of satellite galaxy of Ml 06. 

(2) We consider galaxies with no resolved structure 
but showing rather extended 10^^) and faint surface 
brightness structure (extremely faint surface bright- 
ness) as candidate satellite galaxies. 

(3) Galaxies with similar radial velocity to that of 
Ml 06 are considered to be satellites. 

Finally we find 16 satellite candidates and five 



probable satellite candidates of Ml 06. Hereafter, we 
call the candidate satellite galaxies as satellite galax- 
ies. Twelve of the satellite galaxies are previously 
known, and four of them are new findings. We list the 
properties of these satellite galaxies in Table 3. The 
galaxy type is estimated based on the morphology of a 
galaxy on CCD images and the surface brightness pro- 
files. We display an atlas of 16 satellite galaxies in Fig. 
2. Relatively large satellite galaxies (S2, S3, S4, S7, 
and S15) show elliptical shapes. Two intermediate-size 
galaxies, S5 and SIO are typical edge-on disk galaxies. 
S13 has a structure of typical irregular galaxy, being 
close to a very bright foreground star. We masked out 
the sky regions of this bright star before deriving the 
surface photometry. Two satellite galaxies S6 and S9 
showing spherical shapes with relatively low central 
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Figure 4. (g — r) surface colour profiles of the satellite galaxies in M106. Solid lines represent the box-car smoothed profiles 
of the observed data. Dashed lines represent the linear fits of the surface colour profiles for the intermediate radial range 
(see Table 4 for details), in each panel. The typical error of the surface colour ranges from 0.1 to 0.2. The meaning of the 
vertical arrows and lines are the same as those in Fig. 3. 



surface brightness seem to be dwarf spheroidal galax- 
ies. Surface brightness profiles of satellite galaxies S8 
and S9 with spherical appearance are well represented 
by exponential profiles for the entire radial ranges. 
Two satellite galaxies showing spherical shapes with 
extremely low surface brightness, Sll and S14, are 
morphologically similar to dwarf spheroidal galaxies. 
S16 shows an elongated shape and its surface bright- 
ness profiles for the outer region are well represented 
by an exponential law. S 16 is a dwarf elliptical galaxy 
with a central nucleus. SI was located at the bound- 
ary of CCD chip, so that we could not estimate its 
morphological type. The surface photometry of SI in- 
dicates that the morphological type of SI is close to 
a dwarf elliptical galaxy. We classify PI as probable 
satellite galaxy of Ml 06 because its central surface 
brightness is much fainter than all the satellite galax- 
ies. Remaining galaxies, P2, P3, P4, and P5, which 



have extremely small effective radii (< 3.5 arcsec), are 
also classified as probable satellite galaxies of Ml 06. 

3.2 Surface Brightness and Color Profiles 

We display the radial surface brightness profiles of all 
the satellite galaxies of M106 in Fig. 3, and list the 
results of surface photometry including (g — r) colour 
gradient in Table 4. The radial profiles of surface mag- 
nitude in g and r bands are displayed in open circles 
and open squares, respectively. We marked the posi- 
tion of the effective radius (reff) within which half of 
the total galaxy luminosity is emitted by thin arrows, 
together with the positions of two characteristic radii, 
r25 and rnoim, by the vertical solid lines and the verti- 
cal dotted lines. r25 and rHoim are defined as the radii 
where the B-band surface brightness becomes 25 mag 
arcsec"^ and 26.5 mag arcsec"^, respectively. 
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Table 3. A list of the candidate satellite galaxies in M106 



ID 




R.A. 




Dec. 


^out 


Cout 


PAout 


9^ 




Vr 


Name 


Type 


Selection^ 


SI 


12 


14 


05.00 


47 


26 08.6 


12.9 


0.329 


146 


19.20±0.06 


18.47±0.05 


- 


- 


dE 


2 


S2 


12 


15 


46.41 


48 


07 51.0 


39.8 


0.294 


48 


13.49±0.06 


12.88±0.05 


766^ 


NGC4218 


Sd 


1,3 


S3 


12 


16 


11.71 


47 


52 59.7 


124.5 


0.636 


50 


11.88±0.04 


10.95±0.03 


920^ 


NGC4220 


Sa 


3 


S4 


12 


17 


49.90 


47 


24 33.0 


113.0 


0.553 


20 


12.86±0.24 


12.12±0.16 


484^ 


NGC4248 


Im 


1,3 


S5 


12 


18 


11.07 


46 


55 00.7 


41.0 


0.725 


110 


16.29±0.06 


16.29±0.09 


431^ 


- 


Sm 


1,3 


S6 


12 


19 


06.21 


47 


43 49.8 


10.3 


0.100 


170 


19.30±0.06 


18.74±0.08 


- 


KK132 


dSph 


1 


S7 


12 


19 


09.39 


47 


05 27.5 


48.1 


0.250 


70 


15.37±0.03 


14.75±0.03 


202^ 


UGC7356 


dE,N 


1,3 


S8 


12 


19 


33.21 


47 


27 05.3 


16.7 


0.111 


170 


17.32±0.22 


16.76±0.44 


780^ 


KK134 


dE 


2 


S9 


12 


19 


35.96 


47 


05 35.2 


11.1 


0.000 





18.82±0.10 


18.13±0.05 






dSph 


2 


SIO 


12 


20 


17.49 


48 


08 15.6 


55.6 


0.727 


97 


15.95±0.04 


15.56±0.03 


797^ 


UGC7392 


Sd 


1,3 


Sll 


12 


20 


30.07 


47 


29 24.8 


18.5 


0.450 


140 


20.10±0.54 


19.42±0.45 






dSph/dIrr 


1 


S12 


12 


20 


40.18 


47 


00 03.2 


22.5 


0.231 


80 


17.41±0.03 


16.69±0.02 




KK136 


dE 


2 


S13 


12 


20 


48.42 


47 


49 33.4 


39.3 


0.667 


115 


15.66±0.20 


15.23±0.21 


762^ 


UGC7401 


Im 


1,3 


S14 


12 


20 


54.86 


46 


49 48.8 


13.1 


0.214 


45 


19.22±0.16 


18.49±0.16 






dSph 


2 


S15 


12 


23 


27.94 


46 


59 37.8 


121.1 


0.539 


187 


11.69±0.04 


10.96±0.04 


770^ 


NGC4346 


SO/SBO 


1,3 


S16 


12 


23 


46.18 


47 


39 31.9 


22.3 


0.417 


140 


18.39±0.09 


17.82±0.09 






dE,N 


1 


PI 


12 


18 


18.76 


47 


48 17.0 


6.5 


0.000 





19.86±0.16 


19.86±0.12 






dsph 




P2 


12 


18 


06.54 


46 


28 50.1 


2.7 


0.188 


135 


20.96±0.17 


21.07±0.22 






dIrr/dSph 




P3 


12 


17 


06.33 


46 


54 08.2 


2.2 


0.146 


130 


19.94±0.07 


20.19±0.07 






dlrr/dSph 




P4 


12 


17 


02.25 


46 


54 20.2 


4.1 


0.475 


172 


19.77±0.06 


19.87±0.06 






dlrr/dSph 




P5 


12 


16 


37.40 


47 


09 12.0 


1.8 


0.000 





20.33±0.10 


20.69±0.23 






dIrr/dSph 





"Applied criteria for individual candidate satellite galaxies. 
^An average value of NED and SDSS. 
^Radial velocity from NED. 
^Radial velocity from SDSS. 



Most of the satellite galaxies of Ml 06 show sur- 
face brightness profiles typical of a disk galaxy, while 
the surface brightness profiles of S15 resemble those 
of an elliptical galaxy. Surface brightnesses of three 
satellite galaxies S6, Sll, and S14 are extremely low 
so that we could not determine the two characteristic 
radii, r25, and rnoim for these galaxies. These low sur- 
face brightness galaxies are expected to be an analogy 
of Sextans and Sagittarius dwarf spheroidal galaxies 
in the Local Group. Radial surface brightness profiles 
of the satellite galaxies, S7, S9, and S16 show a strong 
signature of the presence of the central nucleus. The 
radial surface brightness profile of S7 near the galaxy 
centre is found to be very close to that of the point 
source profile, implying the presence of the central nu- 
cleus. On the other hand, the central structures of S9 
and S16 are slightly extended compared with the typ- 
ical point source profile. 

We fit the derived radial surface brightness pro- 
files with a single exponential law (thin solid lines 
over data points) except for S15, for which we use the 
de Vaucouleures law. The physical properties of the 
satellite galaxies of Ml 06 derived by fitting the radial 
surface brightness profiles are listed in Table 4. The 
radial surface brightness profiles of the outer regions 
of the satellite galaxies of Ml 06 are fitted well by an 
exponential profile. However, near the galaxy centre 
there are seen a variety of structures. We divide the 
satellite galaxies of Ml 06 with exponential disk into 
4 categories based on their central surface brightness 



profiles: (1) galaxies with surface brightness profiles 
represented well by a single exponential law (SB and 
S12), (2) galaxies expected to host a clear /marginal 
central nucleus (S7, S9, and S16), (3) satellite galaxies 
showing brighter central surface brightness than the 
exponential law but more extended than the PSF (S3, 
S4, S5, Sll, and S13), and (4) galaxies with fainter 
central brightness than the exponential disk profiles 
(SI, S2, S6, SIO and S14). 

The central structures of the satellite galaxies in 
category 3 might be explained by the presence of the 
central bulge in these galaxies. The appearances of 
these galaxies are, however, quite different from one 
another. Moreover, two galaxies, S5 and SIO, showing 
typical shapes of disk edge have quite different central 
structure. 

We present the derived radial {g — r) colour pro- 
files of the satellite galaxies of M106 in Fig. 4. For 
clear appearance we apply the box car smoothing ker- 
nel to the data (solid lines), {g — r) surface colour 
of the satellite galaxies ranges from 0.1 to 1.0 mag 
arcsec"^, corresponding to 0.2 and 1.2 mag arcsec"^ 
in (B — V). We note that the expected foreground red- 
deni ng toward M1Q6 is negligible {E{B - V^) < 0.02, 
Schl egel et al.l (| 19981 )1. We plot the radial colour pro- 
files of 16 satellite galaxies using the same scale in 
y-axis for clear comparison. We derive the colour gra- 
dient of {g — r) colour profiles of these galaxies and 
summarized the results in Table 4. For the linear fit 
of the surface colour profile we only include the data in 
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Table 4. Surface Photometry of satellite galaxies in M106 



ID 




eit 


R^ 


nolm 


^9 




Fit ranges^ 


{g — r) slope^ 


SI 


23.69^ 


6.3 


3.5 


8.9 


4.09±0.17 


3.99±0.16 




,1 n 


0.018±0.119 


S2 


19.63^ 


11.9 


30.8 


40.6 


6.82±0.10 


7.11±0.11 




-^zU 


-0.035±0.078 


S3 


17.46^ 


32.0 


99.3 


134.6 


22.88±0.40 


22.89±0.40 




,1 nn 
iUU 


-0.147±0.066 


S4 


20. 24^ 


43.9 


86.0 


128.7 


27.02±0.30 


28.55±0.33 






0.057±0.040 


S5 


21. 43^ 


15.3 


21.8 


34.5 


9.54±0.21 


10.05±0.24 






0.176±0.090 


S6 


25.00^ 


8.1 


- 


7.4 


6.85±0.30 


6.42±0.27 




,1 o 


-0.115±0.172 


S7 


20. 71^ 


23.9 


24.8 


44.9 


16.76±0.46 


16.30±0.44 


0.2- 


■^50 


-0.041±0.046 


S8 


23.18^ 


11.6 


7.6 


17.3 


7.25±0.10 


7.27±0.10 


0.2- 


■^40 


0.045±0.053 




oq 7q2 


7 7 


n 7 


7 




c Qo-i-n 1 A 

o.yuzizu. 1^ 


0.2- 


-^25 


— n n9R-i-n n7i 


SIO 


21. 74^ 


15.5 


27.3 


40.4 


9.28±0.19 


9.88±0.22 


0.3- 


■^30 


0.069±0.075 


Sll 


25.59^ 


18.5 




0.7 


16.43±0.63 


18.40±0.80 


0.2- 


■^30 


0.054±0.061 


S12 


23.50^ 


14.4 


6.8 


18.7 


9.19±0.14 


9.11±0.14 


0.2- 


■^30 


0.003±0.061 


S13 


21. 25^ 


30.9 


33.5 


70.2 


22.46±0.87 


20.34±0.69 


2.5- 


■^40 


0.193±0.140 


S14 


25.44^ 


11.4 




7.6 


7.49±0.26 


8.91±0.36 


0.3- 


■^25 


-0.016±0.067 


S15 


16.70 


21.5 


91.6 


135.3 






10.0- 


■^150 


-0.060±0.151 


S16 


22. 12^ 


13.5 


1.8 


14.0 


11.91±0.61 


12.35±0.65 


2.0- 


■^30 


0.041±0.139 


PI 


25.77 


9.2 




1.3 


5.39±0.34 


3.87±0.17 


0.5- 


■^20 


-0.328±0.074 


P2 


23.80 


1.8 




4.9 


1.03±0.04 


0.73±0.02 


0.5^ 


-4 


-0.489±0.239 


P3 


22.99 


2.3 


3.2 


5.1 


1.71±0.04 


1.60±0.03 


0.5- 


■^10 


-0.213±0.111 


P4 


22.80 


2.5 


4.2 


5.8 


1.60±0.03 


1.60±0.03 


0.5- 


■^10 


-0.156±0.111 


P5 


24.59 


3.4 


1.3 


5.3 


2.90±0.24 


2.43±0.17 


0.5^ 


-6 


-0.301±0.170 




denotes the central surface bri^ 


^htness in units of mag arcsec 


-2 







^Effective radius(i^eff )5 standard radius(i?25), and Holmberg radius (i^holm) measured in r-band 
images are given in units of arcsec. 

^Disk scale length(if) in g and r bands in units of arcsec. 

^Radial ranges for linear fit of {g — r) colour profiles in units of arcsec. 

^{9 ~ f") colour gradient in units of mag/Alogr where r in units of arcsec. 

1,2,3,4 Categories for surface brightness profiles and central surface brightness as described in 



the intermediate radial range as summarized in Table 
4. The surface colour profiles of S3 and S15 get bluer 
with increasing radius, while those of S5 and S13 get 
redder as radius increases. The radial surface colour 
profiles of remaining 12 galaxies do not show any clear 
signature of colour gradient. Despite the similar trend 
in the radial surface brightness profiles, the colour pro- 
files of S3 and S5 differ clearly: slow increase in {g — r) 
colour with increasing radius near the central region 
followed by a slower decrease (getting bluer) for S3, 
but the trend is quite opposite for S5. The marginally 
visible resolved stars in S5 might indicate the pres- 
ence of young MS stars and/or AGB stars, which are 
brighter than the tip of the red giant branch stars. 
This indicates the recent star formation in S5. For 
three galaxies (S6, Sll, and S14) with extremely low 
surface brightness, we found no trend in the colour 
profiles due to the large photometric errors. 



3.3 Spatial Distribution 

The two dimensional distribution of the 16 satellite 
galaxies of M106 (see Fig. 1) shows (1) the central 
concentration of the satellite galaxies toward Ml 06 
and (2) the higher probability to find satellite galaxies 
along the major axis of Ml 06. We investigate these 



features in more detail using quantitative analysis as 
shown in Fig. 5. 

We show one dimensional spatial distributions of 
the 16 satellite galaxies in Fig. 5. Fig. 5 (a) and (c) 
present the surface number density profiles as a func- 
tion of the projected distance from the centre of M106 
and from the major axis of Ml 06. Both profiles show 
clear signatures of concentration toward the centre of 
M106 and toward the major axis of M106, respectively. 
We fit the radial number density profile of the satel- 
lite galaxies using a power- law for the radial range of 
^ 45 arcmin, finding a power index of —2.12 ± 0.55. 
A power index of a power- law profile for Fig. 5(c) is 
— 1.28 ± 0.52 much flatter compared to the power in- 
dex for galacto- centric radial distribution (Fig. 5(a)). 
The derived power-law index for the galacto- centric 
surface number density profile is very similar to the 
expected value for isothermal distribution. This im- 
plies that 12 of 16 candidate satellite galaxies with 
dMioG < 50 arcmin 100 kpc) found in this study 
are probably genuine satellite galaxies of Ml 06. 

Because we are dealing with small number statis- 
tics {N < 20), the derived power index of power- law 
for the number density profiles might be affected sig- 
nificantly by the binning method. To examine the ef- 
fect of binning for the number density profiles we de- 
rive the number density profiles using four different 
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Figure 5. Spatial distributions of the satellite galaxies 
of M106: (a) number density distribution against the pro- 
jected distance from the centre of M106, (b) cumulative 
number distribution of the galacto-centric radial distance 
of the satellite galaxies from M106, (c) number density 
distribution of the satellite galaxies against the projected 
distance from the major axis of M106, and (d) cumu- 
lative number distribution with increasing projected dis- 
tance from the major axis of M106. The solid lines in Figs. 
5(a) and 5(c) are the best-fitting power-law for the ranges 
< dMiOG < 45/, and < dMioe-plane < 30/. The shaded 
regions in Figs. 5(b) and 5(d) represent the expected cu- 
mulative number distribution with la boundary assuming 
a spatially uniform distribution. 

binning methods; (1) an equal number of data points 
in each bin, (2) binning the data points in an equal 
linear interval, (3) binning the data points in an equal 
logarithmic interval, and (4) binning the data points 
in an equal r^^^ interval. We applied four different 
binning methods to obtain the radial number den- 
sity profiles of 16 satellite galaxies of Ml 06. We then, 
fit the data with c^m loe < 45 arcmin to the power- 
law for the each density profile, obtaining power in- 
dices of -2.90 ±0.89, -2.14 ±0.74, -1.93 ±0.51, and 
—2.21 ± 0.59, respectively. All of the power indices 
agree well within the errors, and the derived power- 
index (—2.12) is close to the weighted mean of the 
values. 

The satellite galaxies in the outer region (> 100 
kpc) show more or less like a flat distribution rather 
than the power-law. We fit the radial number den- 
sity profile of the satellites including the galaxies in 
the outer region of M106 (c/mio6 > 100 kpc), obtain- 
ing a power index of —1.33 ± 0.36. This is somewhat 
shallower than the power index for the radial number 
density profile of 12 satellite galaxies. However due to 
the large uncertainty it is quite uncertain whether the 




Total V magnitude 



Figure 6. Distribution of both the candidate satellite 
galaxies (open circles) and the probable satellite galaxies 
(star symbols) of M106 in total magnitude versus central 
surface brightness (Eq) in V-band. The solid and dashed 
lines represent the best-fit linear relations for the data of 16 
candidate satellite galaxies and 16+5 satellite galaxies, re- 
spectively. The hashed region around the solid line denotes 
la boundary of the best-fit linear relation. 



satellite distribution of the outer region of M106 is flat 
or not. One of the causes for the flat distribution of 
Ml 06 satellite located at the outer region may be the 
inclusion of the background galaxies in the list of the 
Ml 06 satellite galaxies. 

Fig. 5 (b) and (d) show the cumulative number 
distribution of satellite galaxies against the projected 
distance from the centre of M106, and that of the pro- 
jected distance from the major axis of Ml 06, respec- 
tively (solid lines). For comparison we display the ex- 
pected cumulative number distribution of 16 artificial 
sources assuming spatially uniform distribution with 
shaded regions. The width of the shaded regions rep- 
resents la boundary. The concentrated distribution 
of the satellite galaxies toward the centre of Ml 06 is 
also shown clearly with a high statistical confidence 
(> 95%). For the outer region (dMioe > 50 arcmin) 
the observed distribution of the satellite galaxies is 
very similar to that of the uniform distribution (This 
is also observed in Fig. 5 (a)). In addition to the cen- 
tral concentration toward the centre of Ml 06, they are 
located preferentially along the major axis of Ml 06. 
This implies that a large portion of the satellite galax- 
ies of Ml 06 were formed around Ml 06 and are cap- 
tured by Ml 06 moving along the merging filaments, 
which are spatially aligned with the major axis of 
M106. 
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Figure 7. Cumulative luminosity function of the satellite 
galaxies in M106 (star symbols), LG (open triangles), and 
the M81 group (open circles). The solid lines overlaid on 
the data-points are the best-fitting cumulative Schechter 
functions. The photometric bands for each data-set are 
noted in the figure. 

4 DISCUSSION 

4.1 Misidentification and Missing Satellite 
Galaxies 

The sample of 16 satellite galaxies of Ml 06 found 
in this study might contain non members. Most of 
the misidentifications are expected to be connected 
with faint background galaxies. The probability to in- 
clude background galaxies as the satellite galaxies of 
M106 can be higher in applying the criterion (3) in 
§3.1 because the background galaxies (but still within 
a few Mpc from M106) with large peculiar motion 
could possibly be identified as the satellite galaxies of 
Ml 06. In Table 3, there are 9 galaxies with known ra- 
dial velocities. Considering the radial velocity of M106 
(vr — 448kms~"^), S3 that shows the largest radial 
velocity difference (At'r = 472km~^) still seems to 
be within the expected dynamical boundary of M106 
{\/\vr{sat - host)] ^ 400 - SOOkms"^). Thus we con- 
clude that the probability including misidentification 
in 9 satellite galaxies with radial velocity information 
is negligible. 

For the remaining 7 galaxies without radial ve- 
locity information, we are not able to discriminate 
the genuine satellites from the misidentified ones. We 
need to measusre their radial velocities by spectro- 
scopic observations or to derive their distances by high 
resolution photometry. Observations with higher spa- 
tial resolution using space telescope and/or ground- 
based observations equipped with adapted optics at 



near-infrared wavelengths are crucial to confirm these 
galaxies as the satellite galaxies of Ml 06, by detecting 
a significant number of red giant branch stars in the 
satellite galaxies. 

One of the possible contaminations in the current 
list of satellite galaxies of Ml 06 is the inclusion of 
background galaxies especially for the outer region of 
M106. The number of satellite galaxies at (Imiog > 50 
arcmin is four. Two (S2 and SIO) of these 4 galaxies 
have the radial velocity of Vr ~ 750 km/sec, somewhat 
larger than the radial velocity of Ml 06 but still smaller 
than or similar to the radial velocity of S3 and S8. 
For the remaining two galaxies (SI and S16) without 
radial velocity information we are not able to decide 
if they are dynamically bound to M106. 

To examine the misidentification of background 
galaxies as the satellite galaxies of Ml 06, we show the 
central surface brightness versus total magnitude of 
both candidate satellite galaxies (open circles) and 
probable satellite galaxies (star symbols) in Fig. 6. 
Two candidate satellite galaxies (S4 and S16) are 
located somewhat outside of the la boundary, but 
still within ~ 2a boundary. For the probable satellite 
galaxies, we find that only PI is located well inside 
the boundary. The fitting error for the best-fit linear 
relation using all the data (16+5 galaxies) is ^ 40% 
larger than that for the best-fit linear relation using 16 
candidate satellite galaxies. S16 shows some resolved 
stars with an extended halo, which satisfies the most 
strong selection criterion. S4 (NGC4248) also shows 
some resolved stars, and the radial velocity (484 km 
s~^) is very similar to that of Ml 06. Therefore, we 
conclude that both S4 and S16 are candidate satellite 
galaxies of Ml 06. 

In addition to the misidentification of the satellite 
galaxies of Ml 06, there might be a number of galaxies 
not included in the current list of satellite galaxies of 
Ml 06 because of incomplete photometry. We estimate 
the expected number of missing satellite galaxies due 
to the gaps in the present CFHT mosaic observations 
(see §2 for the description of gaps). The width of the 
each gap is ^ 0.9 arcmin, resulting in 360 arcmin^ in 
the entire observed field. This corresponds to ~ 3% of 
the whole observational field of view. Since the gap is 
rather narrow the probability of missing a large (size 
> 1 — 2 arcmin) satellite galaxy is considerably small. 
We estimate the expected number of missing small 
satellite galaxies (r ^ 1 arcmin) located in these gaps 
as follows: First we assume a uniform distribution of 
the satellite galaxies. We then assume that if a half 
of a galaxy is located in this gap we fail to find this 
galaxy as the satellite galaxy of Ml 06. The expected 
number of missing satellite galaxies due to the gaps 
of mosaic observation is ^ 0.7. Because the distribu- 
tion of the satellite galaxies of Ml 06 is concentrated 
toward the centre of M106 and the distribution of the 
gaps is roughly uniform, the number might be smaller 
than 0.7. We compared the catalogue of dwarf galax- 
ies around M106 of Karachentsev et al. (2007) with 
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Figure 8. (B — V)o versus My diagram of the satellite 
galaxies of M106 (open circles with error bars). As a com- 
parison the data for the Local Group dwarf galaxies are 
displayed with other symbols. The solid line represents a 
boundary to separate spirals or irregulars from dSph or 
Ellipticals (lMateolll998h . 

our list of satellite galaxies. We found that all dwarf 
satellite galaxies of Karachentsev et al. (2007) near 
Ml 06 are common in our list except for the dwarf 
galaxy, dl2 17+4703 that is found to be located in one 
of the gaps. This is consistent with the expected num- 
ber of missing satellite galaxies due to the gaps be- 
tween CCDs in the present study. We do not include 
this satellite galaxy, dl 2 17+4703 for the subsequent 
analysis because we do not have its morphology and 
photometric properties. 

It is more difficult to estimate the number of miss- 
ing satellite galaxies not located in these gaps. One 
of the causes that make it difficult to identify the 
faint satellite galaxies of Ml 06 is the observational 
limits; poor spatial resolution and shallow photomet- 
ric depth. Identification of satellite galaxies of M106 
based on the present observations is mainly dependent 
on surface brightness rather than the total magnitude. 
Two satellite galaxies (Sll and S14) have low central 
surface brightness (/io,v ^ 25 mag arcsec"^). There- 
fore we conclude that satellite galaxies with extremely 
faint surface brightness may be missing in the present 
list of satellite galaxies of Ml 06. This is more severe 
for the galaxies located in F3 because the measured 
sky brightness of F3 is much brighter than Fl and F2. 

We discuss qualitatively the limit of our de- 
tec tion method by c ompa ring our results th those 
of IChiboucas et al.l (|2009l ) because their observa- 
tions were also performed using CFHT MegaCam. 
IChiboucas et al. 1 f2009) found 22 new dwarf galaxies 




Figure 9. Spatial distribution of the satellite galaxies 
of M106 (star symbols), MW (open triangles), M31 (filled 
triangles), and the M81 group (open circles) against (a) 
sky-projected radial distance to the centre of each galaxy, 
and (b) sky-projected distance to the major axis of each 
galaxy. To mimic the data-set of M106 we transform the 
galacto-centric Cartesian coordinates of the satellite galax- 
ies of Milky Way, M31, and the M81 group assuming an 
inclination angle of 67° for M106. 

in the M81 group using a larger field (65 square de- 
gree) MegaCam survey of the M81 group, resulting 
in a total of 44 galaxies in the M81 group includ- 
ing M81. To find the satellite galaxies they chose two 
different methods; visual inspection by two of three 
authors, and a 2-point correlation auto-detection rou- 
tine. They detected 17 out of 22 new candidates of 
the M81 group by visual inspection of the MegaCam 
images. To find the satellite galaxies missed in the 
visual inspection they applied a 2-point correlation 
auto-detection routine for the photometric catalogue 
of point sources. The auto-detection method calcu- 
lates the 2-point correlation function for the point 
sources with 22 < < 27 on scales of 0.02 — 0.11 kpc, 
then a spatial region with a higher c orrelation values 
might be a candidate satellite galaxy. IChiboucas et al.l 
(|2009l ) applied this method for their photometric cat- 
alogue of the M81 group because they were able to de- 
tect a large number of TRGB stars in the dwarf galax- 
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ies of the M81 group. Assuming Mi^trgb — —4.05 
and typical colours of TRGB of {V — I)trgb 1.5, 
and {V - R)trgb - 0.8 (|Le3ll993l l. Mr' rob is de- 
tected at ~ —3.1, corresponding to t'^^rgb ^ 24.7 at 
the distance of M81. They argued that since their ob- 
servational limit {r' ^ 25.0 — 25.5) is slightly fainter 
than the TRGB brightness they detected a significant 
amount of TRGB stars in the satellite galaxies of the 
M81 group. Auto detection method adds 5 more satel- 
lite galaxies to the catalogue of satellite galaxies of the 
M81 group, mostly galaxies with faint surface bright- 
ness. 

To investigate the limit of our detection method 
we compare the observational parameters of both 
MegaCam observations. Our MegaCam observations 
of M106 is rather deeper (exposure times in r-band is 
^ 2 times longer than r^-band observations of the M81 
group). The atmospheric seeings of the both observa- 
tions are similar to ^ 0.75 arcsec, and both observa- 
tions were carried out at dark nights. Therefore, the 
photometric limiting magnitude of the present obser- 
vation is expected to be ^ 0.75 mag fainter than the 
MegaCam observations of the M81 group. However, 
because the distance to Ml 06 is ^ 2 times farther 
than the distance to M81 {d ^ 3.6Mpc), TRGB is ex- 
pected to be detected at r ^ 26.2 for our MegaCam 
survey, similar to or fainter than the magnitude limit 
of the present MegaCam observation of Ml 06. There- 
fore, we expect to detect no or an extremely small 
number of TRGB stars in satellite galaxies of Ml 06. 
Since 2-point correlation auto-detection method needs 
a catalogue of points sources in the satellite galaxies 
we are not able to apply this method for our Mega- 
Cam observations. Our catalogue of satellite galaxies 
of Ml 06 does not include the galaxies with extremely 
faint central surface brightness, which were detected in 
the M81 group (4 out of 22 galaxies with fio{r) > 27). 
We might miss those extremely faint galaxies in the 
present catalogue of satellite galaxies of Ml 06. 



4.2 Comparison with LG and the M81 group 

MW and M31 in LG have a large number of satellites 
whose properties are well known. The M81 group is 
known to have 43 satellite galaxies of which about half 
are new findings (Chiboucas et al. 2009). We compare 
the properties of the satellite galaxies of M106 with 
those of MW, M31, and the the M81 group. Because 
the brightness of 16 satellite galaxies of Ml 06 is much 
brighter than those of the extremely faint galaxies of 
MW, we only consider the brighter {My ^ —9) satel- 
lite galaxies of MW (samples in Mateo 1998; van den 
Bergh 2000; Metz et al. 2007). 

We show the cumulative luminosity function (LF) 
of satellite galaxies of Ml 06 compared to those of 
the satellite galaxies in LG and the M81 group in 
Fig. 7. For M106 and the M81 group we show the 
r-band and r^-band cumulative LFs while we display 
y-band cumulative LF for LG galaxies. We exclude 



the error bars of data points in Fig. 7 for the clear 
comparison. Then we fit the observed cumulative LFs 
to cumulative Schechter LFs (solid lines). To com- 
pare the LFs of satellite galaxies we do not include 
the host galaxies of each system (MW, M31, M81, 
and M106). For the M81 group, we also exclude the 
faintest sample from the fit to the Schechter LF, which 
is significantly fainter than the majority of satellite 
galaxies of the M81 group. The faint end slopes of 
Schechter LFs for M106, LG, and the M81 group 

are -1.19tH6, -l-Oet^s, and -1.29tHL respec- 
tively. The faint end slo pe of the M81 group is basi- 
cally identical to that of IChiboucas et ahl (|2009l l. The 
faint end slope of Schechter LF of M106 is slightly 
steeper than that of LG but flatter by ~ 0.1 com- 
pared to the M81 group. The cumulative LF of the 
M81 group rises suddenly at M^/ — 13 due to the 
inclusion of a large n umber of satellite ga laxies found 
in MegaCam survey (|Chiboucas et al ][^9). A signif- 
icantly larger number of satellite galaxies of the M81 
group compared to Ml 06 might be due to inclusion of 
some a rtifacts and/or b ac kgro und galaxies as pointed 
out bv lChiboucas et al.l (|2009h (See the comments of 
their Table 1). They performed spectroscopic confir- 
mation of only 3 BCD candidates among 22 new satel- 
lite candidates, which turned out to be the members 
of the M81 group. The cumulative LFs of LG and 
Ml 06 satellite galaxies are represented well by the 
cumulative Schechter function, while the cumulative 
Schechter LF does not represent well the observed cu- 
mulative LF of the M81 group as shown in Fig. 7. The 
faint end slope of cumulative LF of Ml 06 is much flat- 
ter than the expected faint end slope by ACDM model 
of -1.8 (Trentham k Tully, 2002). 

In Fig. 8, we display {B — V)o — My colour- 
magnitude diagram of 16 satellite galaxies in Ml 06 
(open circles with error bars). We transform SDSS g 
and r magnitudes to Johnson B and V magnitud es us- 
ing the trans formation relation bv | Chonis Gaskelll 
([2008,). 

Since IChiboucas et al.1 (|2009l ) have observed 
the M81 group using only r^-band, we show only the 
colour-magnitude distribution of M106 and LG satel- 
lites. The distribution in the colour-magnitude dia- 
gram (CMD) for LG dwarf galaxies associated with 
neither MW nor M31 are shown as a comparison. 
The solid line is a boundary based on observations 
to divide the dwarf galaxies into early type galaxies of 
dSph or ellipticals (redder of the solid line) and late 
type galaxies of spirals or irregulars (Mateo 1998). 
The locations of dwarf satellite galaxies of Ml 06 in 
the CMD coincide roughly with early type galaxies of 
M31 and MW. However the CMD distribution of LG 
dwarf galaxies not related with large spiral galaxies is 
quite different from that of Ml 06 satellites. 

The mean colour of dwarf satellite galaxies of 
M106 is estimated to he < [B — V)q >~ 0.73, consis- 
tent with that of the dwarf satellite galaxies of MW 
and M31. Interestingly two brightest satellite galax- 
ies of M106 (S3 and S15) have much redder {B - V) 
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Figure 10. Distribution of dwarf galaxies of M106 and 
LG in disk scale length (H) versus My - For dwarf galaxies 
of LG we include dirr and dSph. The solid line represents 
the relation for the dSph galaxies in LG (|van den Berghl 
[2OQQI 



colour than the dwarf galaxies of MW with similar lu- 
minosity. This difference is related with the morpho- 
logical difference between them. The brightest satel- 
lite galaxy of Ml 06 is an early type spiral galaxy, while 
that of MW is a late type spiral. James et al.J (|2008h 
reported that it is quite rare to find star forming late 
type spiral galaxies as a satellite galaxy of a giant spi- 
ral galaxy like MW. However, there is a high probabil- 
ity of morphology confo rmity between host and satel- 
lite galaxies ( Ann et a l] l2008l ) and similar conformity 
is present between a target galaxy and its companion, 
especially when the companion is close enough to be 
located insi de th e virial radius of the target galaxy 
(jPark et al.ll20Q8h . 

The spatial distribution of the satellite galaxies in 
M106 shows two distinct features; (1) galacto- centric 
concentration toward M106 and (2) preferred concen- 
tration toward the major axis of Ml 06. We compare 
the spatial distribution of the satellites of M106 to 
those of the dwarf satellite galaxies in LG and the M81 
group and display the results in Fig 9. The galacto- 
centric Cartesian coordinates of the satellite galaxies 
of Milky Way and M31 are from Metz et al. (2007) 
and Koch & Grebel (2006), respectively. To mimic the 
sky-projected two dimensional distribution of M106, 
we transform their Cartesian coordinates by rotating 
along the semi-minor axis of the disk. We determine 
the inclination angle of Ml 06 using the known major 
and minor-axis radii assuming zero-thickness of the 
disk of M106 following Ma et al. (1997), obtaining 
~ 67°. Due to the limited survey area of the present 



study, the distribution of the satellite galaxies of M106 
is truncated at ^ 130kpc, corresponding to ^ 1°. 

The surface number density profiles against the 
sky-projected galacto- centric distance (Fig. 9(a)) are 
very similar to one another among four satellite sys- 
tems with an indication of slightly higher number den- 
sity of the satellite galaxies of the M81 group com- 
pared to MW, M31, and M106. The surface number 
density of the M81 group is higher by a factor of 
~ 2.5 than that of LG satellite galaxies for for the 
radial ranges of 50 ^ 200kpc. We fit the radial sur- 
face number density profiles of the satellites of MW, 
M31, and the M81 group with the power- law, obtain- 
ing the power indices of -1.29 =b 0.14, -1.36 ± 0.17, 
and —1.62 zb 0.12, respectively. For the brighter sam- 
ple of the M81 group with M^/ < —10, we obtain 
the power index of —1.24 db 0.18. The radial surface 
number density profiles of the satellites of MW, M31, 
and the M81 group decrease gradually roughly follow- 
ing the power-law with varying power indices. On the 
contrary, the radial distribution of satellites of Ml 06 
shows some hint of flat distribution for the outer re- 
gion (c/miog lOOkpc). M106 is expected to be lo- 
cated in rather denser environment than MW, M31, 
and M81, this might affect the distribution of satel- 
lite galaxies of M106 compared to those of MW, M31. 
and the M81 group. However, since the number den- 
sity profile of Ml 06 satellites is only available for the 
radial ranges of dMioe < 130kpc, we are not able to 
conclude that if the outer distribution of M106 satel- 
lites follows the uniform distribution. To reconcile this 
issue, it is need to extend the survey area to cover 
the whole dynamical boundary of Ml 06. We show 
the surface number density of the satellite galaxies of 
M106, MW, M31, and the M81 group with increasing 
projected distance from the semi-major axis of each 
galaxy in Fig. 9(b). The distributions of the satellite 
galaxies of these galaxies are similar within the un- 
certainties. We also compute the number of satellite 
galaxies of Ml 06 by assuming that the flat distribu- 
tion in Fig. 5(a) (c/mio6 > 45 arcmin) is mainly due 
to background galaxies. We find that the number of 
satellite galaxies located between 10 — 45 arcmin from 
M 106 is ^ 8 by assuming flat background distribution. 
This number might be the lower limit of the actual 
satellite galaxies of Ml 06, corresponding to some half 
of candidate satellite galaxies found in this study. 

To estimate the number of satellite galaxies of 
M106 located outside of the surveyed region (c/miog > 
130kpc) of Ml 06, we count the number of satellite 
galaxies outward of this distance of MW, M31, and 
the M81 group by assuming that the virial mass of 
each system is similar to each other. We find the num- 
bers 4(- 30%), 8(- 50%), and 18(- 40%) for MW, 
M31, and the M81 group, respectively, with an aver- 
age of 40% of total satellite galaxies in these galaxies. 
This simply implies that there might be ^ 10 satel- 
lite galaxies with Mr < —10 in the outer region of 
Ml 06. This estimation might be an upper limit be- 
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cause the virial mass of Ml 06 might be somewhat 
smaller than MW, M31, and the M81 group and the 
catalogue of the satellite galaxies of the M81 group 
includes some very faint galaxies down to M^/ ~ — 7. 
We compute the expected number of satellite galaxies 
of Ml 06 located beyond the surveyed region of Ml 06 
(<^Mio6 > 130kpc) by extrapolating the measured 
power-law slope of —2.12, resulting in N(> 130kpc) 
^ 6. This number is quite similar to that of M31, but 
significantly lower than those of MW 11), and the 
M81 group(^ 21). Thus, to find the expected number 
of satellite galaxies in Ml 06 we need to increase the 
survey area by a factor of ^ 7. 

We show the distribution of dwarf galaxies of 
M106 in Mv versus disk scale length {H) diagram 
in Fig. 10 (open circles). For comparison the distri- 
bution of LG dwarf galaxies are shown with differ- 
ent symbols. We do not include the d ata for the M81 
group because IChiboucas et al.] (|2009l ) presented only 
the effective radii of the satellite galaxies. The scale 
length data for LG galaxies including spirals, dirr, and 
dSph are from van den Bergh (2000). Most of dwarf 
galaxies of Ml 06 are scattered in the domain of the 
distribution of dSph of LG. Sll classified as dSph in 
this study has relatively large disk scale length than 
the other dSph galaxies with similar luminosity. Late- 
type spiral satellite galaxies of M106 (S2, S3 and S4) 
are seen to lie well above the relation for dSph galax- 
ies. Overall the distribution of dwarf galaxies in Ml 06 
is very similar to that of dIrr and dSph of LG in this 
diagram. 



5 SUMMARY AND CONCLUSION 

We present the result of a wide-area survey for satel- 
lite galaxies around M106. Our survey is based on the 
large field deep g and r images obtained using the 
CFHT MegaCam so that we could perform a system- 
atic searching for satellite galaxies of Ml 06. We sum- 
marize the main conclusions of the present study as 
follows. 

1. Based on 1.7° x 2° deep observations of Ml 06, we 
found 16 candidate satellite galaxies of Ml 06, located 
as far as ^ 130kpc from the centre of M106. Four 
of them are new findings. Due to the limited expo- 
sure times and the large atmospheric seeing we are 
able to detect only a small number of TRGB stars in 
the satellite galaxies of Ml 06. Our survey area might 
be significantly smaller than the physical boundary 
of Ml 06 dark halo. Therefore we might miss several 
satellite galaxies outside of ^ 130kpc from M106. 

2. A satellite system of Ml 06 shows a clear signature 
of central concentration toward Ml 06 and is more con- 
centrated toward the semi-major axis of Ml 06. We fit 
the inner region of the radial surface number density 
profiles with a power-law, resulting a power index of 
—2.1 zb 0.5, similar to that of the isothermal distribu- 
tion. 

3. The surface brightnesses of 15 satellites galaxies 



are found to fit well with an exponential profile while 
the remaining brightest satellite galaxy of M106 (SI 5) 
shows a radial profile of the typical early type galaxy. 

4. The luminosity of 16 satellite galaxies spans a large 
range (AMy ^ 7.6 mag), but {B — V) colours have a 
mean value of 0.73 with a small dispersion. The lumi- 
nosity function of the satellite galaxies of Ml 06 is well 
represented by the Schechter function with a faint end 
slope of —1.19. 

5. The spatial distribution of the satellite galaxies of 
M106, Milky Way, M31, and the M81 group is found 
to be similar to one another. 
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